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Abstract 

Double-differential cross sections for light-ion (p, d, t, 3 He and a) production 
in oxygen induced by 96 MeV neutrons are reported. Energy spectra are mea- 
sured at eight laboratory angles from 20° to 160° in steps of 20°. Procedures 
for data taking and data reduction are presented. Deduced energy-differential 
and production cross sections are reported. Experimental cross sections are 
compared to theoretical reaction model calculations and experimental data 
at lower neutron energies in the literature. The measured proton data agree 
reasonably well with the results of the model calculations, whereas the agree- 
ment for the other particles is less convincing. The measured production 
cross sections for protons, deuterons, tritons and alpha particles support the 
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trends suggested by data at lower energies. 

PACS numbers: 24.10.-i, 25.40.-h, 25.40.Hs, 25.40.Kv, 28.20.-v 

1 Introduction 

Fast-nucleon induced reactions provide useful means to investigate nuclear 
structure, to characterize reaction mechanisms and to impose stringent con- 
straints on nuclear model calculations. Although oxygen is a light nucleus 
with doubly closed shells, it can be expected that many statistical assump- 
tions hold for nucleon-induced reactions at several tens of MeV, because the 
level density at high excitation energies is sufficiently high that shell effects 
and other nuclear structure signatures are washed out. Light nuclei also 
have low Coulomb barrier, implying that the suppression of charged-particle 
emission is weak. Therefore, nuclear reaction models for equilibrium and 
pre-equilibrium decay can be tested and benchmarked. Experimental data 
on reactions in oxygen in the literature at incident neutron energies of 27, 
40, and 60 MeV [U |2] and between 25 and 65 MeV S U offer possibilities 
to test the predictions of reaction models. 

In recent years, an increasing number of applications involving fast neu- 
trons have been developed or are under consideration, e.g., radiation treat- 
ment of cancer [3 13 El, neutron dosimetry at commercial aircraft alti- 
tudes jSj, soft-error effects in computer memories [TOl ITT], accelerator-driven 
transmutation of nuclear waste and energy production ^2j EI] , and determi- 
nation of the response of neutron detectors ^1] . Data on light-ion production 
in light nuclei such as carbon, nitrogen and oxygen are particularly important 
in calculations of dose distributions in human tissue for radiation therapy at 
neutron beams, and for dosimetry of high energy neutrons produced by high- 
energy cosmic radiation interacting with nuclei (nitrogen and oxygen) in the 
atmosphere [HI • When studying neutron dose effects in radiation therapy 
and at high altitude, it is especially important to consider oxygen, because 
it is the dominant element (65 % by weight) in average human tissue. 

In this paper, we present experimental double-differential cross sections 
(inclusive yields) for protons, deuterons, tritons, 3 He and alpha particles pro- 
duced by 96 MeV neutrons incident on oxygen. The measurements have been 
performed at the cyclotron of The Svedberg Laboratory (TSL), Uppsala, us- 
ing the MEDLEY experimental setup JH]- Spectra have been measured at 
8 laboratory angles, ranging from 20° to 160° in 20° steps. Extrapolation 
procedures are used to obtain coverage of the full angular range, and con- 
sequently energy-differential and production cross sections are deduced, the 
latter by integrating over energy and angle. The experimental data are com- 
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pared to results of calculations with nuclear reaction codes and to existing 
experimental data at lower incident neutron energies. 

The experimental methods are briefly discussed in Sec. |2] and data re- 
duction and correction procedures are presented in Sec. 01 The theoretical 
framework is summarized in Sec. |U In Sec. EJ the experimental results are 
reported and compared with theoretical and previous experimental data. 
Conclusions and an outlook are given in Sec. |H1 

2 Experimental methods 

The experimental setup, procedures for data reduction and corrections have 
been described in detail recently |T71 ITSj and therefore only brief summaries 
are given here. 

The neutron beam facility at The Svedberg Laboratory (TSL) uses the 
7 Li(p,n) 7 Be reaction to produce a quasi-monoenergetic neutron beam [TH] . 
The lithium target was 8 mm thick in the present experiment and enriched 
to 99.98 % in 7 Li. The 98.5±0.3 MeV protons from the cyclotron impinge on 
the lithium target, producing neutrons with a full-energy peak of 95.6±0.5 
MeV with a width of 1.6 MeV (FWHM). With a beam intensity of 5 //A, the 
neutron flux in the full-energy peak is about 5-10 4 neutrons/(s-cm 2 ) at the 
target location. The collimated neutron beam has a diameter of 80 mm at the 
location of the target, where it is monitored by a thin film breakdown counter 
(TFBC) [20] ■ Relative monitoring was obtained by charge integration of the 
proton beam in a Faraday cup located in the proton beam dump. The two 
beam monitor readings were in agreement during the measurements. 

The charged particles are detected by the MEDLEY setup [16 . It consists 
of eight three-element telescopes mounted inside a 100 cm diameter evacuated 
reaction chamber. Each telescope consists of two fully depleted AE silicon 
surface barrier detectors and a CsI(Tl) crystal. The thickness of the first 
AE detector (AEi) is either 50 or 60 fim, while the second one (AE2) is 
either 400 or 500 /im. They are all 23.9 mm in diameter (nominal). The 
cylindrical CsI(Tl) crystal, 50 mm long and 40 mm in diameter, serves as the 
E detector. 

A 22 mm diameter 500 /zm thick (cylindrical) disk of Si02 is used as the 
oxygen target. For the subtraction of the silicon contribution, measurements 
using a silicon wafer having a 32-32 mm 2 quadratic shape and a thickness of 
303 /im are performed. 

For absolute cross section normalization, a 25 mm diameter and 1.0 mm 
thick polyethylene (CH 2 ) n target is used. The np cross section at 20° labora- 
tory angle provides the reference cross section [21]. Instrumental background 



3 



is measured by removing the target from the neutron beam. It is dominated 
by protons produced by neutron beam interactions with the beam tube and 
reaction chamber material, especially at the entrance and exit of the reaction 
chamber and in the telescope housings. Therefore, the telescopes at 20° and 
160° are most affected. 

The time-of-flight (TOF) obtained from the radio frequency of the cy- 
clotron (stop signal for TDCs) and the timing signal from each of the eight 
telescopes (start signal) is registered for each charged-particle event. Typical 
count rates for target-in and target-out runs were 10 and 2 Hz, respectively. 
The dead time of the data acquisition system was typically 1 — 2 % and never 
exceeded 10 %. 

3 Data reduction procedures and corrections 

The AE — E technique is used to identify light charged particles ranging 
from protons to lithium ions. Good separation of all particles is obtained 
over their entire energy range and particle identification is straightforward. 

Energy calibration of all detectors is obtained from the data themselves [TTf 
ITH] . Events in the AE — E bands are fitted with respect to the energy de- 
posited in the two silicon detectors. This energy is determined from the 
detector thicknesses and tabulated energy loss values in silicon [22] • The 
AEi detectors are further calibrated and checked using a 5.48 MeV alpha 
source. For the energy calibration of the CsI(Tl) detectors, two parameteri- 
zations of the light output versus energy of the detected particle [IB1 HZl HBj 
are used, one for hydrogen isotopes and another one for helium isotopes. Sup- 
plementary calibration points are provided by the H(n,p) reaction, as well 
as transitions to the ground state and low-lying states in the 12 C(n,d) 11 B, 
16 0(n,d) 15 N and 28 Si(n,d) 27 Al reactions. The energy of each particle is ob- 
tained by adding the energy deposited in each element of the telescope. 

Low-energy charged particles are stopped in the AE% detector leading to 
a low-energy cutoff for particle identification of about 3 MeV for hydrogen 
isotopes and about 8 MeV for helium isotopes. The helium isotopes stopped 
in the AEi detector are nevertheless analyzed, and a remarkably low cutoff, 
about 4 MeV, can be achieved for the experimental alpha-particle spectra. 
These alpha-particle events could obviously not be separated from 3 He events 
in the same energy region, but the yield of 3 He is about a factor of 30 smaller 
than the alpha-particle yield in the region of 8 MeV, where the particle iden- 
tification works properly. The assumption that the relative yield of 3 He is 
small is supported by the theoretical calculations in the evaporation peak re- 
gion. In conclusion, the 3 He yield is within the statistical uncertainties of the 
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alpha-particle yield for alpha energies between 4 and 8 MeV. A consequence 
of this procedure is that the 3 He spectra have a low-energy cutoff of about 8 
MeV. 

Knowing the energy calibration and the flight distances, the flight time 
for each charged particle from target to detector can be calculated and sub- 
tracted from the measured total TOF. The resulting neutron TOF is used 
for selection of charged-particle events induced by neutrons in the main peak 
of the incident neutron spectrum. 

Background events, collected in target-out runs and analyzed in the same 
way as target-in events, are subtracted from the corresponding target-in runs, 
with SiC>2 and silicon targets, after normalization to the same neutron flu- 
ence. 

Due to the finite target thickness, corrections for energy loss and parti- 
cle loss are applied to both targets individually. Details of the correction 
methods are described in Refs. EHJ- The cross sections for oxygen are 
obtained after subtraction of the silicon data from the SiC>2 data with proper 
normalization with respect to the number of silicon nuclei in the two targets. 1 

Even if a great majority of the neutrons appears in the narrow full-energy 
peak at 95.6 MeV, a significant fraction (about 13 %) belongs to a tail ex- 
tending towards lower energies, remaining after the TOF cut, see Fig.[0 The 
average neutron energy with the tail neutrons included is 94.0 MeV. The par- 
ticle spectra have not been unfolded with the neutron energy distribution, 
because it is anticipated that the energy variation of the cross sections is 
rather weak in the energy range of interest. Furthermore, the data set is 
called 96 MeV (95.6 MeV) data, because the peak of the distribution is quite 
dominant and any structure observed at the high-energy end of the ejectile 
spectra is due to the peak of the neutron energy dsitribution. The np cross 
section is, however, measured at the peak of the distribution (95.6 MeV) and 
corrected for the tail contribution. The correction to 94.0 MeV is performed 
using the known energy dependence of the np cross section. 

Other corrections of the data are performed in analogy with the similar 
experiment dealing with silicon and described in detail in the corresponding 
publication J7j. The data and method for the efficiency correction of the 
CsI(Tl) detectors, reported in Ref. [TH] and used in Ref. [T7] and the present 
work, have recently [2S] been corroborated by Monte Carlo calculations. 

Absolute double-differential cross sections are obtained by normalizing 
the oxygen data to the number of recoil protons emerging from the CH 2 

1 In the process of extracting the oxygen data, the silicon data of Ref. ^2] have been 
reanalyzed. In doing so, we have adapted some changes and also found two mistakes. See 
Ref. [22 i n this issue. 
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target. After selection of events in the main neutron peak and proper sub- 
traction of the target-out and 12 C(n,p) background contributions, the latter 
taken from a previous experiment, the cross section can be determined from 
the recoil proton peak, using np scattering data [21]. All data have been 
normalized using the np scattering peak in the 20° telescope. 

4 Theoretical models 

The present data have been compared with nuclear theory predictions, com- 
puted with the two nuclear reaction codes GNASH [221 121] and TALYS [2H] • 
While GNASH has been widely used during the last years, TALYS is a new 
code that has just been released in the public domain. The GNASH calcula- 
tion is performed at a neutron energy of 100 MeV with parameters given in 
a recent evaluation for medical purposes j2H] as described in Ref. [12] . Since 
oxygen is at the boundary of the mass range the TALYS code is aimed for, 
the code is described in some detail below. 

Both GNASH and TALYS integrate direct, pre-equilibrium, and statis- 
tical nuclear reaction models into one calculation scheme and thereby give 
predictions for all the open reaction channels. Both codes use the Hauser- 
Feshbach model for sequential equilibrium decay and the exciton model for 
pre-equilibrium emission, though GNASH uses the one-component model, 
i.e. without isospin distinction of the excited nucleons, and TALYS uses the 
two-component model, see below. The angular distributions are obtained 
using the Kalbach systematics jSHj. 

The purpose of TALYS is to simulate nuclear reactions that involve neu- 
trons, photons, protons, deuterons, tritons, 3 He and alpha particles in the 1 
keV - 200 MeV energy range. Predicted quantities include integrated, single- 
and double-differential cross sections, for both the continuum and discrete 
states, residue production and fission cross sections, gamma-ray production 
cross sections, etc. For the present work, single- and double-differential cross 
sections are of interest. To predict these, a calculation scheme is invoked 
which consists of a direct + pre-equilibrium reaction calculation followed by 
subsequent compound nucleus decay of all possible residual nuclides calcu- 
lated by means of the Hauser-Feshbach model. 

For the optical model potentials (OMP) of both neutrons and protons on 
16 O up to 200 MeV, the global OMP of Ref. [SI] was used. These potentials 
provide the necessary transmission coefficients for the statistical model calcu- 
lations. Although the global neutron OMP has been validated for A > 24, at 
the high incident energy considered in this work, an adequate description of 
the basic scattering observables is expected, at least for the incident neutron 
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channel and the high energy inelastic scattering and charge-exchange leading 
to discrete states and the continuum. For the low-energy outgoing charged 
particles, the non-validated use of the global OMP may have larger conse- 
quences. Obviously, a system of a total of 17 nucleons can hardly be called 
statistical, and this short-coming may be reflected in the prediction of some 
of the observables that concern low emission energies. For complex particles, 
the optical potentials were directly derived from the nucleon potentials using 
the folding approach of Watanabe [22] • Finally, since applying the charged- 
particle OMP's for nuclides as light as 16 O may be physically dubious, we 
renormalize the obtained OMP transmission coefficients with the empirical 
non-elastic cross sections of Ref. 

The high-energy end of the ejectile spectra are described by pre-equili- 
brium emission, which takes place after the first stage of the reaction but 
long before statistical equilibrium of the compound nucleus is attained. It is 
imagined that the incident particle step-by-step creates more complex states 
in the compound system and gradually loses its memory of the initial energy 
and direction. The default pre-equilibrium model of TALYS is the two- 
component exciton model |3*4*1 I35j. A remark similar to that given above for 
the OMP applies: the two-component exciton model for nucleon reactions has 
been tested, rather successfully, against basically all available experimental 
nucleon spectra for A > 24 [34J. The current system A = 17, falls outside that 
mass range, and does not entirely qualify as a system that can be handled 
by fully statistical models such as the exciton model. 

We recall the basic formula of Ref. jSl] for the exciton model cross section, 

eg eq 
Ptt Pv 

Pv, h 

where Ptt{Pu) is the proton (neutron) particle number and ^(h^) the proton 
(neutron) hole number, a is the compound formation cross section, and 
S pre is the time-integrated strength which determines how long the system 
remains in a certain exciton configuration. The initial proton and neutron 
particle numbers are denoted p° = Z p and p° u = N p with Z P (N P ) being the 
proton (neutron) number of the projectile. In general, h n = p^ — p® and 
h u = Pv — P„, so that the initial hole numbers are zero, i.e. h® = h° u = 0, 
for primary pre-equilibrium emission. The pre-equilibrium part is calculated 
by Eq. (0), using p^ q = p eq = 6, whereas the remainder of the reaction flux 
is distributed through the Hauser-Feshbach model. In addition, the never- 
come-back approximation is adopted. 
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The emission rate w k for ejectile k with spin s k is given by 
, , i Tp \ 2s k + 1 tp u(p n - Z k ,h n ,p u - N k ,h u ,E x ) 

(2) 

where (Jk,inv(Ek) is the inverse reaction cross section as calculated from the 
optical model, and u is the two-component particle-hole state density. The 
full reaction dynamics that leads to Eq. (JTJ) is described in Refs. j^H • 
We here restrict ourselves to the formulae given above since they contain the 
model- and parameter-dependent quantities. The expression for S pre con- 
tains the adjustable transition matrix element M 2 for each possible transi- 
tion between neutron-proton exciton configurations. A proton-neutron ratio 
of 1.6 for the squared internal transition matrix elements was adopted to give 
the best overall agreement with experiment, i.e., M 2 U = M 2 n = 1.6M 2 n = 
1.6M 2 U = 1.6M 2 . For 16 0, we use the following expression for the matrix 
element 



„, 2 °-6 
M = —r 

A 3 



4.2 x 10 5 

6.8 + 



(3) 



;^ + io.7) 3 _ 

where n is the exciton number. Partial level density parameters g„ = Z/17 
and g u = N/17 were used in the equidistant spacing model for the partial 
level densities. Finally, an effective surface interaction well depth V = 12 
MeV fgj was used. 

At incident energies above several tens of MeV, the residual nuclides 
formed after binary emission may have so large excitation energy that the 
presence of additional fast particles inside the nucleus becomes possible. The 
latter can be imagined as strongly excited particle-hole pairs resulting from 
the first binary interaction with the projectile. The residual system is then 
clearly non-equilibrated and the excited particle that is high in the contin- 
uum may, in addition to the first emitted particle, also be emitted on a short 
time scale. This so-called multiple pre-equilibrium emission forms an alter- 
native theoretical picture of the intra-nuclear cascade process, whereby the 
exact location and momentum of the particles are not followed, but instead 
the total energy of the system and the number of particle-hole excitations 
(exciton number). In actual calculations, the particle- hole configuration of 
the residual nucleus after emission of the ejectile, is re-entered as initial con- 
dition in Eq. When looping over all possible residual configurations, the 
multiple pre-equilibrium contribution is obtained. In TALYS, multiple pre- 
equilibrium emission is followed up to arbitrary order, though for 96 MeV 
only secondary pre-equilibrium emission is significant. 

It is well-known that semi-classical models, such as the exciton model, 
have always had some problems to describe angular distributions (essentially 
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because it is based on a compound-like concept instead of a direct one). 
Therefore, as mentioned previously, the double-differential cross sections are 
obtained from the calculated energy spectra using the Kalbach systemat- 
ics (SOI- 

To account for the evaporation peaks in the charged-particle spectra, 
multiple compound emission was treated with the Hauser-Feshbach model. 
In this scheme, all reaction chains are followed until all emission channels are 
closed. The Ignatyuk model jHHj has been adopted for the total level density 
to account for the damping of shell effects at high excitation energies. 

For pre-equilibrium reactions involving deuterons, tritons, 3 He and alpha 
particles, a statistical contribution from the exciton model is automatically 
calculated with the formalism described above. It is, however, well known 
that for nuclear reactions involving projectiles and ejectiles with different 
particle numbers, mechanisms like stripping, pick-up and knock-out play an 
important role and these direct-like reactions to the continuum are not cov- 
ered by the exciton model. Therefore, Kalbach has developed a phenomeno- 
logical contribution for these mechanisms 37], which is included in TALYS. 
Among the advantages over the older method (which is included in GNASH) 
we mention here a better consideration of the available phase space through 
normalized particle-hole state densities and a better empirical determination 
of the pick-up, stripping, knock-out strength parameters, enabled by the 
more extensive experimental database that is now available. It has recently 
been shown (see Table I of Ref. [38] ) that for medium and heavy nuclides this 
method gives a considerable improvement over the older methods. The latter 
seemed to consistently underpredict neutron-induced reaction cross sections 
involving pick-up of one or a few nucleons. In this paper, the two methods 
meet again, this time for the prediction of reactions on a light nucleus, and 
their performance will be compared in the next section. 

5 Results and discussion 
5.1 Experimental results 

Double-differential cross sections of 16 0(n,xlcp) reactions, where lcp stands 
for light charged particle, at laboratory angles of 20°, 40°, 100° and 140° for 
protons, deuterons, tritons, 3 He and alpha particles are shown in Figs. |2]— 
respectively, all angles plotted with the same cross section scale for each 
emitted particle to facilitate comparison of magnitudes. The choice of energy 
bin width depends on the energy resolution in the experiment, the thick 
target correction and acceptable statistics in each energy bin. The error bars 
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in Figs. El— El represent statistical uncertainties only. 

The overall relative statistical uncertainties of individual points in the 
double-differential energy spectra at 20° are typically 8 % for protons, 13 
% for deuterons, 20 % for tritons, 15 % for 3 He and 12 % for alpha parti- 
cles. As the angular distributions are forward-peaked, these values increase 
with angle. The systematic uncertainty contributions are due to thick target 
correction (1 — 20 %), collimator solid angle (1 — 5 %), beam monitoring 
(2 — 3 %), number of oxygen nuclei (0.1 %), CsI(Tl) intrinsic efficiency (1 
%), particle identification (1 %) and dead time (<0.1 %). The uncertainty 
in the absolute cross section is about 5 %, which is due to uncertainties 
in np scattering angle, contribution from the low-energy continuum of the 
7 Li(p,n) spectrum to the np scattering proton peak (3 %), reference np cross 
sections (2 %) [2T], statistics in the np scattering proton peak (2 %), carbon 
contribution (0.1 %) and number of hydrogen nuclei (0.1 %). 

From Figs. El— EUit is obvious that the charged-particle emission at forward 
angles from 96 MeV neutron irradiation of oxygen is dominated by proton, 
deuteron and alpha particle channels. The yield of deuterons is about a factor 
of 3 lower than for protons and the spectra of the two other particle types 
studied in this work (tritons and 3 He) are more than an order of magnitude 
weaker. All the spectra have more or less pronounced peaks at low energies 
(below 10 — 15 MeV), the angular distributions of which are not too far 
from isotropy except for alpha particles, where the yield at backward angles 
is about four times weaker than at 20°. The low-energy peak is not fully 
observed in the 3 He spectra due to the 8 MeV low-energy cutoff discussed in 
Sec. El 

All the particle spectra at forward angles show relatively large yields at 
medium-to-high energies. The emission of high-energy particles is strongly 
forward-peaked and hardly visible in the backward hemisphere. It is a sign 
of particle emission before statistical equilibrium has been reached in the 
reaction process. In addition to this broad distribution of emitted particles, 
the deuteron spectra at forward angles show narrow peaks corresponding to 
transitions to the ground state and low-lying states in the final nucleus, 15 N. 
These transitions are most likely due to pick-up of weakly bound protons in 
the target nucleus, 16 0. A similar but less pronounced effect is observed in 
the proton spectra at forward angles. The structure observed in this case 
is due to transitions to Gamow- Teller states and other low-lying states with 
considerable single-particle strength pQ. 
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5.2 Comparison with theoretical model calculations 

In Figs. El— IHlthe experimental results are presented together with theoretical 
model calculations. The GNASH calculations of Ref. [2E] have been done 
for protons, deuterons and alpha particles, whereas the TALYS calculations 
discussed in Sec. H] have been performed for all five particle types. The 
TALYS calculations include a transformation of the calculated cross sections 
to the lab system. Also in the GNASH code a similar transformation from 
the cm. to the lab system is performed using the kinematics of one-particle 
emission. Differences between data given in the lab and cm. systems are 
particularly significant in this case, because oxygen is such a light nucleus. 

Fig. |21 shows the comparison between the double-differential (n,px) ex- 
perimental spectra and the calculations based on the TALYS and GNASH 
models. For protons above 25 MeV, both calculations give a reasonably good 
description of the spectra, although the calculated 20° cross sections, in par- 
ticular the TALYS ones, fall below the experimental data. The low-energy 
statistical peak below 15 MeV in the spectra is considerably overpredicted 
by the two codes. The overestimate is particularly strong at backward angles 
for TALYS and at forward angles for GNASH. 

The situation is quite different for the deuteron spectra (Fig. EJ). None 
of the calculations do account very well for the data, although the GNASH 
code gives a reasonable description of the angular dependence of the cross 
section. For the TALYS code deviations between data and calculations of a 
factor of two or more are present. At forward angles the high-energy part is 
strongly overestimated, in particular by the TALYS code, indicating prob- 
lems in the hole-strength treatment. It is obvious, however, that efforts have 
been spent in these calculations to include individual hole-state strengths. 
Such strengths are not included in the GNASH calculations, but in spite 
of this the average behavior of the cross section at high energies is in fair 
agreement with the data. Like for the proton spectra, the statistical peak 
is overpredicted by the TALYS calculations essentially at all angles, whereas 
the GNASH calculations seem to do a slightly better job in this case. 

For tritons (Fig.HJ), the TALYS calculation gives a fairly good description 
of the experimental data, except that it fails to account for an intensity bump 
around 15 MeV observed at forward angles. 

The general trends of the forward-angle 3 He data (Fig. 03) are reasonably 
well described in the TALYS calculations although the cross sections are 
underestimated by a large factor. At backward angles the yield is very small 
and it is difficult to make quantitative comparisons. 

The overall shapes of the alpha particle spectra (Fig. ^ are reasonably 
well described by the two models. The GNASH calculations, however, over- 
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predict the cross sections at forward angles and underpredict them at large 
angles, whereas the TALYS calculations do the opposite, i.e., underpredict 
at small angles and overpredict at large angles. 

The ability of the models to account for the low-energy peak caused by 
evaporation processes (and for a particles also 3a breakup of 12 C) is not 
impressive. In general, the models tend to overpredict the cross sections. It 
should, however, be kept in mind that the peak maximum is close to (for 
3 He below) the low-energy cutoff, which complicates the comparison. An- 
other complication in this context is that the GNASH cross sections although 
given in the lab system, are calculated using the kinematics of one-particle 
emission |2l)l |2*7| for the c.m.-to-lab transformation, which obviously is an 
approximation. 

Experimental angular distributions at low, medium and high ejectile en- 
ergies are shown in Figs. 171— ITT1 for protons, deuterons, tritons, 3 He and alpha 
particles, respectively. The angular distributions are fitted by a simple two 
parameter function, aexp(6cos#) [3D]. The data are compared with angular 
distributions calculated on the basis of the GNASH and TALYS models. In 
general, the TALYS model gives a weaker angular dependence than the data, 
whereas the GNASH model, although being closer to the data, tends to give 
a slightly steeper angular variation. 

A conspicuous deviation from the experimental angular distribution is 
seen for the TALYS prediction at the lowest outgoing energies, e. g. at 
8 — 12 MeV in Fig. [7J We think this is attributed to wrong partial spectrum 
contributions to the total spectrum. The slightly forward-peaked angular 
distribution suggests that the spectrum at these emission energies is not as 
compound-dominated as the TALYS calculation suggests. Instead secondary, 
and even tertiary, pre-equilibrium emission may not be negligible even in the 
evaporation peak. Multiple pre-equilibrium emission is taken into account 
in TALYS but only contributes at somewhat higher emission energies. A 
way to make multiple pre-equilibrium (processes) relatively more important 
is to reduce the compound nucleus emission contribution, but we find that 
the predicted evaporation peak is rather insensitive to parameter variations. 
Hence, this is an open problem for TALYS, which apparently has been solved 
for the GNASH calculation. 

5.3 Integrated spectra 

For each energy bin of the light-ion spectra, the experimental angular dis- 
tribution is fitted by a simple two-parameter function, aexp(6cos#) [3U], as 
exemplified in the previous subsection (Figs. UJ— fTTjl . This allows extrapola- 
tion of double-differential cross sections to very forward and very backward 
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angles. In this way coverage of the full angular range is obtained. By inte- 
gration of the angular distribution, energy-differential cross sections (da/dE) 
are obtained for each ejectile. These are shown in Fig. El together with the- 
oretical calculations. For all ejectiles both calculations give a fair description 
of the energy dependence. Both calculations are in good agreement with the 
proton experimental data over the whole energy range, although the calcu- 
lation for (n,p) reactions to discrete states is underestimating the data. A 
study of the spectroscopic strengths for these states would be welcome. Con- 
cerning the deuteron spectra, the GNASH calculations are in good agreement 
with the data, whereas the TALYS code gives cross sections a factor of two 
or more larger than the experimental ones at energies above 30 MeV. In 
the case of alpha particles, the GNASH calculation tends to overpredict the 
high-energy part of the spectrum, and the TALYS calculations fall below the 
data above an alpha particle energy of 25 MeV. The energy dependence of 
the triton and 3 He spectra are well described by the TALYS code, but in 
both cases the calculation falls below the data above about 20 MeV. 

The production cross sections are deduced by integration of the energy- 
differential spectra (see Table P). To be compared with the calculated cross 
sections, the experimental values in Table Q have to be corrected for the 
undetected particles below the low-energy cutoff. This is particularly im- 
portant for 3 He because of the high cutoff energy. The corrections obtained 
with TALYS seem to be too small in some cases, in particular for the (n,xa) 
production cross section. This is illustrated in Fig. bottom panel, where 
the TALYS curve falls well below the experimental da/dE data in the 4-7 
MeV region. 

The proton, deuteron, triton, and alpha particle production cross sections 
are compared with previous data at lower energies [5] in Fig. ED There 
seems to be general agreement between the trends of the previous data and 
the present data points. The curves in this figure are based on a GNASH 
calculation |29j . 

6 Conclusions and outlook 

In the present paper, we report an experimental data set on light-ion produc- 
tion in oxygen induced by 96 MeV neutrons. Experimental double-differential 
cross sections (d 2 a/dQdE) are measured at eight angles between 20° and 
160°. Energy-differential (da/dE) and production cross sections are obtained 
for the five types of outgoing particles. Theoretical calculations based on nu- 
clear reaction codes including direct, pre-equilibrium and statistical models 
give generally a good account of the magnitude of the experimental cross sec- 
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tions. For proton emission, the shape of the spectra for the double-differential 
and energy-differential cross sections are well described. The calculated and 
the experimental alpha-particle spectra are also in fair agreement with the 
exception of the high energy part, where the GNASH model predicts higher 
yield and the TALYS model lower yield than experimentally observed. For 
the proton evaporation peak, the global TALYS calculation overestimates 
the data. A future activity should be an adjustment of the responsible OMP 
and level density parameters (as was done in the case of GNASH) instead of 
relying on a full global prediction. For the other complex ejectiles (deuteron, 
triton and 3 He) there are important differences between theory and experi- 
ment in what concerns the shape of the spectra at various angles. We think 
this is due to the use of statistical models such as the Hauser-Feshbach model 
and the pre-equilibrium exciton model in mass ranges where these models 
become suspect, and the absence of a break-up model in the theoretical anal- 
ysis. Apart from the aforementioned break-up model, predictions of emission 
of alpha particles may be particularly sensitive to a correct knock-out model 
and the use of adequate complex particle optical model potentials. Stripping 
and knock-out models, level densities, optical model and omission of break- 
up reactions may all add up to problems for something as light as oxygen. 
This needs to be studied in much more detail. Finally, the magnitude of the 
angle- integrated cross sections is reasonably well accounted for. 

For the further development of the field, data at even higher energies 
are requested. The results suggest that the MEDLEY facility, which was 
used in the present work, should be upgraded to work also at 180 MeV, i.e., 
the maximum energy of the TSL neutron beam facility. At present, a new 
neutron beam facility is under commissioning at TSL 39 , covering the same 
energy range, but with a projected intensity increase of a factor five. This 
will facilitate measurements at energies higher than in the present work. 
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Figure captions 



1. Neutron energy distribution with the time-of- flight criterion applied 
derived from np scattering data at an angle of 20 degrees. The peak 
(95.6 MeV), median (95.1 MeV) and average (94.0 MeV) are indicated 
by solid, dashed and dotted vertical lines, respectively. 

2. Experimental double-differential cross sections (filled circles) of the 
0(n,px) reaction at 96 MeV at four laboratory angles. The curves in- 
dicate theoretical calculations based on GNASH (dotted) and TALYS 
(solid). 

3. Experimental double-differential cross sections (filled circles) of the 
0(n,dx) reaction at 96 MeV at four laboratory angles. The curves in- 
dicate theoretical calculations based on GNASH (dotted) and TALYS 
(solid). 

4. Experimental double-differential cross sections (filled circles) of the 
0(n,tx) reaction at 96 MeV at four laboratory angles. The curve indi- 
cates theoretical calculations based on TALYS. 

5. Experimental double-differential cross sections (filled circles) of the 
0(n, 3 Hex) reaction at 96 MeV at four laboratory angles. The curve 
indicates theoretical calculations based on TALYS. 

6. Experimental double-differential cross sections (filled circles) of the 
0(n,ctx) reaction at 96 MeV at four laboratory angles. The curves in- 
dicate theoretical calculations based on GNASH (dotted) and TALYS 
(solid). 

7. Angular distributions of 0(n,px) cross section at ejectile energies of 8 
- 12 MeV (filled circles), 40 - 44 MeV (filled triangles) and 68 - 72 
MeV (open squares). The dashed curves are fits to the data and the 
dotted and solid curves represent calculations based on the GNASH 
and TALYS models, respectively. 

8. Same as Fig. [7| but for the 0(n,dx) cross section. 

9. Same as Fig. [7|but for the 0(n,tx) cross section. No calculations based 
on the GNASH model are available for tritons. 

10. Angular distributions of 0(n, 3 Hex) cross section at ejectile energies of 
10-15 MeV (filled circles), 40 - 45 MeV (filled triangles) and 65 - 70 



19 



MeV (open squares). The dashed curves are fits to the data and the 
solid curves represent calculations based on the TALYS model. 

11. Same as Fig. [TD|but for the 0(n,ax) cross section. The dotted curves 
represent calculations based on the GNASH model. 

12. Experimental energy-differential cross sections (filled circles) for neutron- 
induced p, d, t, 3 He and a production at 96 MeV. The curves indicate 
theoretical calculations based on GNASH (dotted) and TALYS (solid). 

13. Neutron-induced a) proton, b) deuteron, c) triton, and d) alpha particle 
production cross section as a function of neutron energy. The full circles 
are from the present work, whereas the open circles are from previous 
work [Sj. The curves are based on a GNASH calculation. The data 
as well as the calculations correspond to cutoff energies of 6 MeV for 
protons and deuterons and 12 MeV for tritons and alpha particles. 
Note that the cutoff energies are different from those in Table ^ 
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Table 1: Experimental production cross sections for protons, deuterons, tri- 
tons, 3 He and alpha particles from the present work. The experimental data 
in the second column have been obtained with cutoff energies of 2.5, 3.0, 3.5, 
8.0 and 4.0 MeV for p, d, t, 3 He and alpha particles, respectively. The third 
and forth columns show data corrected for these cutoffs, using the GNASH 
(Ref. j2H|) and the TALYS calculations of the present work, respectively. 
Theoretical values resulting from GNASH and TALYS calculations are in 
the fifth and sixth columns, respectively. 

u vro d Experiment Experiment [cutoff corr.] Theoretical calculation 



(mb) GNASH TALYS GNASH TALYS 



(n,px) 


224±11 


248 


231 


259.9 


221.7 


(n,dx) 


72±4 


80 


73 


73.4 


131.3 


(n,tx) 
(n, 3 Hex) 


20±1 




20 




10.6 


6.9±0.6 




8.7 




8.2 


(n,ax) 


132±7 


218 


132 


224.7 


88.4 
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